Microcanonical finite-size scaling in specific heat diverging 2 order phase transitions 
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A Microcanonical Finite Site Ansatz in terms of quantities measurable in a Finite Lattice allows 
to extend phenomenological renormalization (the so called quotients method) to the microcanonical 
ensemble. The Ansatz is tested numerically in two models where the canonical specific-heat diverges 
at criticality, thus implying Fisher-renormalization of the critical exponents: the 3D ferromagnetic 
Ising model and the 2D four-states Potts model (where large logarithmic corrections are known 
to occur in the canonical ensemble). A recently proposed microcanonical cluster method allows to 
simulate systems as large as L — 1024 (Potts) or L — 128 (Ising). The quotients method provides 
extremely accurate determinations of the anomalous dimension and of the (Fisher-renormalized) 
thermal v exponent. While in the Ising model the numerical agreement with our theoretical expec- 
tations is impressive, in the Potts case we need to carefully incorporate logarithmic corrections to 
the microcanonical Ansatz in order to rationalize our data. 

PACS numbers: 05.50.+q. 



I. INTRODUCTION 

The canonical ensemble enjoys a predominant position 
in Theoretical Physics due to its many technical advan- 
tages (convex effective potential on finite systems, easily 
derived Fluctuation-Dissipation theorems, etc.) [33]. This 
somehow arbitrary choice of ensemble is justified by the 
Ensemble Equivalence property, that holds in the Ther- 
modynamic Limit for systems with short range interac- 
tions. 

However, in spite of this long standing prejudice in fa- 
vor of the canonical ensemble, the canonical analysis of 
phase transitions is not simpler. The advantages of mi- 
crocanonical analysis of first-order phase transitions has 
long been known [1, 2], and indeed become overwhelming 
in the study of disordered systems [3]. Furthermore, the 
current interest in mesoscopic or even nanoscopic sys- 
tems, where Ensemble Equivalence docs not hold, pro- 
vides ample motivation to study other statistical ensem- 
bles and, in particular, the microcanonical one [4]. Be- 
sides, microcanonical Monte Carlo [5] is now as sim- 
ple and efficient as its canonical counterpart (even mi- 
crocanonical cluster algorithms are known [2]). Under 
such circumstances, it is of major interest the extension 
to the microcanonical framework of Finite-Size Scaling 
(FSS) [6-9] for systems undergoing a continuous phase 
transition. 

The relation between the microcanonical and the 
canonical critical behavior is well understood only in the 
Thermodynamic Limit. A global constraint modifies the 
critical exponents, but only if the specific-heat of the un- 
constrained system diverges with a positive critical expo- 
nent a > [10] (however, see [11]). The modification in 
the critical exponents, named Fisher renormalization, is 
very simple. Let L be the system size, and consider an 



observable O (for instance, the correlation length) whose 
scaling behavior in the infinite-volume canonical system 
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Now, let e be the internal energy density and e c = 
( e )i=™T ■ Consider the microcanonical expectation 
value of the same observable O in (1), but now at fixed 
energy e. The scaling behavior (1) translates to [34] 
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We will denote the microcanonical exponents with the 
subindcx "m" . Hence, the Fisher renormalization of the 
correlation length exponent v, is v — ► v m = — a), 
that of the order parameter exponent is /3 — > /3 m = (3/(1— 
a), etc. On the other hand, the anomalous dimension rj = 
•q m is invariant under Fisher renormalization [10]. See 
also [12], for a recent extension of Fisher renormalization 
to the case of logarithmic scaling corrections. 

As for systems of finite size, the microcanonical 
FSS [13-15] is at the level of an Ansatz. This Ansatz 
is obtained from the canonical one merely by replacing 
the free-energy density by the entropy density, and using 
Fisher renormalized critical exponents. The microcanon- 
ical Ansatz reproduces the canonical one [16], and it has 
been subject of some numerical testing [15, 17]. Further- 
more, systems undergoing Fisher Renormalization (due 
to a global constraint other than the energy) do seem to 
obey FSS as well [18]. 

A difficulty lies in the fact that the current forms of 
the microcanonical FSS Ansatz [13-15] are in a some- 
how old-fashioned form. Indeed, they are formulated in 
terms of quantities such as e c or the critical exponents, 
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which are not accessible in the absence of an analyti- 
cal solution. In this respect, a great step forward was 
achieved in a canonical context [19] when it was realized 
that the Finite-Lattice correlation length [20] allows to 
formulate the FSS Ansatz in terms of quantities com- 
putable in a Finite-Lattice. This formulation made prac- 
tical to extend Nightingale's phenomcnological renormal- 
ization [21] to space dimensions D > 2 (the so-called 
quotients method [22]). 

Here, we will extend the microcanonical FSS Ansatz to 
a modern form, allowing us to use the quotients method. 
We will test numerically this extended Ansatz in two 
models with a > 0, hence undergoing nontrivial Fisher 
Renormalization, namely the D = 3 ferromagnetic Ising 
model, and the D = 2 four-states ferromagnetic Potts 
model. The Potts model has the added interest of suffer- 
ing, in its canonical form, quite strong logarithmic cor- 
rections to scaling that arc nevertheless under relatively 
strong analytical control [23]. It will be, therefore, quite 
a challenge to control the logarithmic corrections in the 
microcanonical setting. 

The layout of the rest of this paper is as follows. In 
Sect. II, we briefly recall the particular microcanonical 
ensemble used in this work (Lustig's microcanonical set 
up [5] , where the Fluctuation-Dissipation formalism of [2] 
applies). In Sect. Ill we present our extended micro- 
canonical FSS Ansatz. A brief description of simulated 
models and measured observable is presented in Sect. IV 
while the specific simulation details are given in Sect. V. 
The results both for the D = 3 Ising model and for the 
D = 2 Potts model are given in Sect. VI and VII respec- 
tively. Finally we devote Sect. VIII to the conclusions. In 
addition, in Appendix A we propose an extension of the 
quotients method, aimed to speed up convergence to the 
large L limit in the presence of multiplicative logarithmic 
corrections. 



II. THE MICROCANONICAL ENSEMBLE 

The first step in the construction of the ensemble is an 
extension of the configuration space. We add N(= L D ) 
real momenta, pi, to our TV original variables, Oi (named 
spins here) [2, 5]. Note that this extended configuration, 
{o~i,Pi}, appears in many numerical schemes (consider, 
for instance, Hybrid Monte Carlo [24] simulations in Lat- 
tice Gauge Theory) . We shall work in the microcanonical 
ensemble for the {<Ti,??i} system. 

Let U be the original spin Hamiltonian (e.g. Eq. (35) 
in our case). Our total energy is [35] 



N 



£ = J2^T +U (e = £/N , u = U/N). (3) 



The momenta contribution, 
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is necessarily positive, and it is best thought of as a "ki- 
netic" energy. In this mechanical analog, the original spin 
Hamiltonian U can be regarded as a "potential" energy. 
The canonical partition function is (/3 = 1 /T) 
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where X){<r } denotes summation over spin configura- 
tions. Hence, the {pi} play the role of a Gaussian ther- 
mostat. The {pi} are statistically uncorrelated with the 

\ ca: 



spins. Since («;)^ onical = 1/(2/3), one has (e) ca 



\ canonical 



V(2/3). 



Furthermore, given the statistical independence of k 
and u, the canonical probability distribution function for 

e, Pg (e), is merely the convolution of the distributions 
for k and u: 



P { s L) (e) 
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{n)P^ u {e-n). (6) 



In particular, note that for spin systems on a finite lattice, 
Pp L ^' u (u) is a sum of (order TV) Dirac's 8 functions. Now, 
since the canonical variance of k is l/(/3v2TVj, roughly 
s/N discrete u-levels, with u ~ e— 1/(2/3), give the most 

significant contribution to P^ L \e). We see that the mo- 
menta's kinetic energy provide a natural smoothing of 
the comb-like Pi L ^ u (u). Once wc have a conveniently 

smoothed Pi L \e), we may proceed to the definition of 
the entropy. 

In a microcanonical setting, the crucial role is played 
by the entropy density, s(e, TV), given by 



exp[TVs(e,TV)] 
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Integrating out the {pi} using Dirac's delta function in 
(7) we get 



exp[TVs(e, TV)] 



(27rTV)f 
AT(TV/2) 



^>(e, U ,TV), (8) 



uj(e,u,N) = (e — u) ^~9{e — u) 



(9) 



The step function, die — u), enforces e > u. Eq. (8) 
suggests to define the microcanonical average at fixed e 
of any function of e and the spins, 0(e, {o~i}), as [5] 



_ 0(e,{q t }M e , M ,TV) 

We use Eq. (8) to compute ds/de: [2] 
ds(e,TV) 



de 

/3(e;M) 
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Keeping in mind the crucial role of the generating- 
functional in Field-Theory (see e.g. [9]), we extend the 
definition (7) by considering a linear coupling between 
the spins and a site dependent source field hf. 

exp[Ns(e,{hi},N)] = J^Y, ^ ^(iVe-f) , 

(13) 

where £ = Ne is still given by Eq. (3), without including 
the source term. In this way, the microcanonical spin cor- 
relation functions follow from derivatives of s(e, {hi}, N): 
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In particular, if the source term is uniform hi = h we 
observe that the microcanonical susceptibility is given 
by standard fluctuation-dissipation relations, see Rcf. 9 
and Eq. (42), below. 



A. Ensemble equivalence 

Eq. (7) ensures that the canonical probability density 
function for e is 



N 



^' Z N {(3) 
hence, Eq. (11), 



cxp[iV( S ( e ,iV)-/3e)], (15) 



logP^ea) - = N I" de (0) e - 

lei 
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The relation between the canonical and the micro- 
canonical spin- values is given by 

/■CO 

(G) canonical = / fa (Q) e Pjf\e) . (17) 



the saddle point approximation becomes exact: 



\ canonical 
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and we have Ensemble Equivalence: 
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It follows that l/[d0) e /de] at ef^ will tend in the 
large-L limit to minus the canonical specific heat. Thus, 
if the critical exponent a is positive, Eq. (19) will fail 
precisely at e c . Hence, Eq. (21) can be expected to hold 
for all e but e c (or for all f3 but /3 C ). 



B. Double peaked histogram 

The situation can be slightly more complicated if 
Pp c (e) presented two local maxima, remindful of phase 
coexistence. This is actually the case for one of our mod- 
els, the D = 2, four states Potts model [25] . From Eq. (16) 
it is clear that the solution to the saddle point equa- 
tion (18) will no longer be unique. We borrow the fol- 
lowing definitions from the analysis of first-order phase 
transitions (where true phase coexistence takes place) [2] : 

• The rightmost root of (18), e\ is a local max- 
imum of Pp corresponding to the "disordered 
phase" . 

• The leftmost root of (18), e° L a, is a local maximum 
of Pp L>> corresponding to the "ordered phase" . 

• The second rightmost root of (18), e* L * is a local 



minimum 



of P, 



(L) 



Maxwell's construction yields the finite-system critical 
point, f3 CiL (see Fig. 7): 



de 



(22) 



Now, Eqs. (15) and (17) imply that the canonical mean- 
value will be dominated by a saddle-point at e SP , 



(18) 



which can be read as yet another expression of Thermo- 
dynamics second-law, Tds = de . 

The condition of thermodynamic stability (namely 
that (p) e be a monotonically decreasing function of e) 
ensures that the saddle point is unique and that e SP is a 
maximum of P@(e). Under the thermodynamic stability 
condition and if, in the large L limit, 



d(P)e 



de 



< 0. 



(19) 
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and the finite-system estimator of the "surface tension" 
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de (0) e - /3 c ,i 



(23) 



Of course, in the large-L limit and for a continuous tran- 
sition, S L -> 0, (3% -> (3 C and e A L ^ c L , e° L ,0 c , L e c ■ 



III. OUR MICROCANONICAL FINITE-SIZE 
SCALING ANSATZ 

Usually, the Microcanonical FSS Ansatz takes the form 
of a scaling form for the entropy density [13-15]. In 
close analogy with the canonical case, one assumes that 
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s(e, {hg}, N) can be divided in a regular part, and a sin- 
gular term s s ; ng (e, {hg}, N). The regular part is sup- 
posed to converge for large L (recall that N = L D ) to 
a smooth function of its arguments. Hence, all critical 
behavior comes from s s i ng (e, {hg}, N). Note as well that 
we write {hg}, instead of {hi}, to emphasize the spa- 
tial dependence of the sources (supposedly very mild [9]). 
Hence, 

s sing (e,{hg},N) = L- D g(L^{e-e c ),{Ly h hg}) . 

(24) 

Here, g is a very smooth function of its arguments, while 
ijh = 1 + is the canonical exponent, see e.g. [9], 
which does not get Fisher-renormalized. Scaling correc- 
tions due to irrelevant scaling fields, have been ignored 
by other authors [13-15], but will be important for our 
precision tests. We will propose here alternative forms 
of the Ansatz (24), more suitable for a numerical work 
where neither e c nor the critical exponents are known 
beforehand. 

Our first building block is the infinite-system micro- 
canonical correlation length, £oo, e . Indeed, Ensemble 
Equivalence implies that, in an infinite system, the long- 
distance behavior of the microcanonical spin-spin prop- 
agator G(f;e) = {<Jgcrg + f) e - (<Tg) e (ag + i) e behaves for 
large r as in the canonical ensemble (close to a critical 
point £oo,e is large, so that rotational invariance is recov- 
ered in our lattice systems): 



G(f;e) 



.4 



,-r/to 



(25) 



where A is a constant. In particular, note that Ensemble- 
Equivalence implies that the anomalous dimension rj does 
not get Fisher-renormalized. We expect £oo,e = £^- n ° mcal 
if the correspondence between e and T are fixed through 

„ / ~\ canonical 

e — \ e /L=oc,T ■ 

The basic assumption underlying the FSS Ansatz is 
that the approach to the L — > oo limit is governed by the 
dimensionless ratio L/^^. Hence, our first form of the 
Ansatz for the observable O whose critical behavior was 
discussed in Eq. (2) is 



(O) 



L,e 



L—fo(L/^ e ) + 



(26) 



In the above, the dots stand for scaling-corrections, while 
the function fo is expected to be very smooth (i.e. differ- 
cntiable to a large degree or even analytical). A second 
form of the Ansatz is obtained by substituting the scaling 
behavior £00^ oc |e — e c \~ UlTl : 



(Oh 



L 



XJ ^fo(LV»™(e-e c )) + 



(27) 



Again, fo is expected to be an extremely smooth function 
of its argument [36]. In particular, this is the form of 
the Ansatz that follows from Eq. (24) by derivating with 
respect to e or from the source terms. 

However, the most useful form of the Ansatz is ob- 
tained by applying (26) to the Finite-Lattice correlation 



length £l.c, obtained in a standard way (see Ref. [9]) from 
the finite-lattice microcanonical propagator. We expect 
^L,e/L to be a smooth, one-to-one function of L/^^, 
that can be inverted to yield £/£oo,e as a function of 
^L,e/L. Hence, our preferred form of the FSS Ansatz is 
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L~ U G 
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L 



(28) 

Here, Fo and Go are smooth functions of their argu- 
ments and lu is the first Universal Scaling Corrections 
exponent. 

It is important to note that exponent w does not get 
Fisher-renormalized. Indeed, let us consider an observ- 
able O with critical exponent xo at a temperature T 



such 



\ canonical 
'L=oo,T 



Now, ensemble equivalence tells us 



that 0^™,T al = L=oo , e and that ££ a *° ni< 



=oo,T 



Eliminating T in favor of 



O canonical t x O v r a i d £~ L 
L=oo,T = ^L=oc.el A + B Ot, L= 



canonical 



see e.g. [9], we have 
+ ...], (29) 



where A$ and Bq are scaling amplitudes. It follows that 
u> m = uj, and that xq/v — xo,m/v m . 



A. The quotients method 

Once we have Eq. (28) in our hands, it is straightfor- 
ward to generalize the quotients method [22]. In Ap- 
pendix A we describe how it should be modified in the 
presence of (multiplicative) logarithmic corrections to 
scaling. 

Let us compare data obtained at the same value of e 
for a pair of lattices L\ = L and L2 = sL with s > 
1. We expect that a single e Cy L 1 ,L 2 exists such that the 
correlation-length in units of the lattice size coincides for 
both systems: 



sL 



(30) 



Hence, if we compare now in the two lattices the observ- 
able O in (28), precisely at e c ^L,sL, we have 



(O) 



sL.e c 



(0)L, ec , Ll , 



s ^ [1 + A , S L- 



(31) 



where Ao, s is a non-universal scaling amplitude. One 
considers this equation for fixed s (typically s = 2), and 
uses it to extrapolate to L = 00 the L-dependent es- 
timate of the critical exponents ratio xo.m/^m ■ At the 
purely numerical level, mind as well that there are strong 
statistical correlations between the quotients in (30) and 
in (31), that can be used via a jackknife method (see 
e.g. [9]) to strongly reduce the statistical errors in the 
estimate of critical exponents. 

In this work, we shall compute the critical exponents 
from the following operators (x is the susceptibility, while 
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£ is the correlation length, see Sect. IV for definitions): 
X -> x = i/ m (2 - 77) , (32) 

<9 e £ -> 2:0 = + 1 • (33) 

As for the L dependence of e c .L. s , it follows from 
Eq. (27) as applied to £l/L for the two lattice sizes L 
and sL [6, 9]: 

e c , L , s = e c + B \ T S + (34) 

(B is again a non- universal scaling amplitude) . In partic- 
ular, if one works at fixed s, e c ^L,sL tends to e c for large 
L as L- {u>+ ^ ] [37]. 

IV. MODELS AND OBSERVABLES 

We will define here the Model and Observablcs of a 
generic D-dimensional Q-stdXes Potts model. The nu- 
merical study has been done for two instances of this 
model: the three dimensional Ising (Q = 2) model, and 
the two dimensional Q = 4 Potts model. 

We place the spins Uj = 1, . . . , Q at the nodes of a 
hypercubic D-dimensional lattice with linear size L and 
periodic boundary conditions. 

The Hamiltonian is 

u = - E <Wi > ( 35 ) 

<hj> 

where < i, j > denotes first nearest neighbors. For a 
given spin, a, we define the normalized Q- vector s, whose 
q-th component is 




A Q components order parameter for the ferromagnetic 
transition is 

i 

where i runs over all the lattice sites. We will now con- 
sider microcanonical averages. The spatial correlation 
function is 

C(r' - r) = (s(r) ■ s(r')) 

Our definition for the correlation length at a given in- 
ternal energy density e, is computed from the Fourier 
transform of C 



at zero and minimal (||fe m m|| = 2tt/L) momentum [9, 20]: 

V /C(0)/C(fc min ) - 1 
^ L) = 2^jL) ■ (40) 

Note that C can be easily computed in terms of the 
Fourier transform of the spin field, s(fc), as 

C(k)=L D (s(k)-s(-k)) e , (41) 

and that the microcanonical magnetic susceptibility is 

X = L D (M 2 ) e = C(0). (42) 

For the specific case of the Ising model, the traditional 
definitions, using Si = ±1 (recall that Sj = ±l/y/2) are 
related with those of the general model through: 

Rising = _ J2 S.Sj =2U- 3L D , 

<i,3> 

Rising = (43) 

X Isins - 2*. 

Notice that in D = 2 this model undergoes a phase 
transition in (3 C = log(l + y/Q) which is second order for 
Q < 4 and first order for Q > 4 [26]. 

V. SIMULATION DETAILS 

We have simulated systems of several sizes in a suit- 
able range of energies (see Table I). To update the spins 
we used a Swendsen-Wang (SW) version of the micro- 
canonical cluster method [2]. This algorithm depends 
on a tunable parameter, k, which should be as close as 
possible to (/3) e in order to maximize the acceptance of 
the SW attempt (SWA). This requires a start-up using 
a much slower Metropolis algorithm for determining k. 
In practice, we performed cycles consisting of 2 x 10 3 
Metropolis steps, k refreshing, 2 x 10 3 SWA, and a new 
k refreshing. We require an acceptance exceeding 60% 
to finish these pre-thermalization cycles fixing k for the 
following main simulation where only the cluster method 
is used. 

In both studied cases, we have observed a very small 
autocorrelation time for all energy values at every lattice 
size. In the largest lattice for the four states Potts model 
we have also consider different starting configurations: 
hot, cold and mixed (strips). Although the autocorrela- 
tion time is much smaller, for safety we decided to discard 
the first 10% of the Monte Carlo history using the last 
90% for taking measurements. 

VI. RESULTS FOR THE D = 3 ISING MODEL 



(7(fcH]TC(r)c ifc - , 

r 



(39) 



In Fig. 1 (upper panel) we show a scaling plot of 
the correlation length (in lattice size units) against (e — 
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Model 


L 


iV m (xl0 6 ) 


iVo 


Energy range 


Q = 2, D = 3 


8 


20 


42 


[-0.8,-0.9] 




12 


20 


42 


[-0.8,-0.9] 




16 


20 


49 


[-0.8, -0.9] 




24 


20 


25 


[-0.845, -0.875] 




32 


20 


16 


[-0.87, -0.860625] 




48 


20 


10 


[-0.87, -0.860625] 




64 


.-) 


10 


[-0.870625, -0.865] 




96 


5 


10 


[-0.870625, -0.865] 




128 


5 


7 


-0.869375, -0.865625] 


Q = 4, D = 2 


32 


1024 


61 


[-1.2,-0.9] 




64 


128 


61 


[-1.2,-0.9] 




128 


32 


41 


[-1.08,-0.98] 




256 


32 


24 


[-1.08,-1.005] 




512 


25.6 


32 


[-1.07,-1.01] 




1024 


6.4 


30 


[-1.06,-1.02] 



TABLE I: Simulation details for the two considered models. 
For each lattice size L we show the number of measurements 
iV m at each energy and the total number of simulated energies 
uniformly distributed in the displayed energy range N c . For 
the Q = 4, D = 2 model, the value of Nm reported have been 
reached only at specific energies near the peaks of the Maxwell 
construction. Also, additional non-uniformly distributed en- 
ergy values have been simulated near the peaks. 



e c )i 1// " m . For the susceptibility we plot \ ~ L 2 ~ 71 (lower 
panel). If data followed the expected asymptotic criti- 
cal behavior with microcanonical critical exponents they 
should collapse in a single curve. In Fig. 1 we have used 
the canonical critical quantities from Rcfs. 27, 28 trans- 
formed to the microcanonical counterparts using Eq. (2). 
From the plot it is clear that important scaling correc- 
tions exist in both cases for the smallest lattices although 
they are mainly eliminated in the biggest systems. 

To obtain the microcanonical critical exponents we 
used the quotients method, see Sect. Ill A. The clear 
crossing points of the correlation length for different lat- 
tice sizes can be seen in Fig. 2. The determination of 
the different quantities at the crossings, and the position 
of the crossing itself, requires to interpolate the data be- 
tween consecutive simulated energies (see Fig. 2). We 
have found that the method of choice, given de high num- 
ber of energy values available, is to fit, using the least 
squares method, a selected number of points near the 
crossing to a polynomial of appropriate degree. Straight 
lines do not provide good enough fits, however, second 
and third order polynomials give compatible results. In 
practice, we have fitted a second order polynomial using 
the nine nearest points to the crossing, also comparing 
the results with those using the seven nearest points that 
turn out fully compatible. For error determination we 
have always used a jackknife procedure. 

The numerical estimates for e c , £,L,e c /L and the critical 
exponents v m and r], obtained using the quotients method 
for pair of lattices (L, 2L) are quoted in Table II. Our 



^ 0.5 



0.4 



1.2 



0.8 



8 ♦ 48 •''•^^/ > . 

12 >< 64 - '™3<H<* w 

16 * 96 • 

24 128 ♦ + + * 

32 o H 1 1 1- 



♦ 



4 



-l 



-0.5 





(e - e c ) L 



0.5 



FIG. 1: (color online). Scaling plot of the correlation length 
(in lattice size units) and the scaled susceptibility for the 
three-dimensional Ising model. We used the critical values, 
e c = —0.867433 and v m — 0.7077. Notice the strong scaling 
corrections for the small systems, as well as the data collapse 
for the largest lattices. 



L 


e c ,L,2L 








8 


-0.861831(12) 


0.44922(3) 


0.8033(42) 


0.0564(2) 


12 


-0.865010(10) 


0.46106(5) 


0.7968(31) 


0.0492(4) 


16 


-0.866020(6) 


0.46710(5) 


0.7717(22) 


0.0469(4) 


24 


-0.866767(3) 


0.47411(4) 


0.7665(11) 


0.0437(3) 


32 


-0.867034(4) 


0.47813(6) 


0.7594(13) 


0.0425(5) 


48 


-0.867228(2) 


0.48278(5) 


0.7492(5) 


0.0412(3) 


64 


-0.867302(2) 


0.48555(11) 


0.7457(16) 


0.0397(8) 



TABLE II: Lattice size dependent estimates of critical quan- 
tities for the microcanonical D = 3 Ising model. The dis- 
played quantities are: crossing points e c ,i,2£ for the corre- 
lation length in units of the lattice size, £/L itself at those 
crossing points, and the estimates for the correlation length 
exponent v m and the anomalous dimension n. All quantities 
are obtained using parabolic interpolations. 



small statistical errors allow to detect a tiny L evolution. 
An extrapolation to infinite volume is clearly needed. 

Before going on, let us recall our expectations as ob- 
tained applying Fisher renormalization to the most accu- 
rate determination of canonical critical exponents known 
to us [i/ m = v/{l — a) = vj{Dv — 1)]: 

u m = 0.7077(5) (from v = 0.6301(4) [29]) , (44) 
??m = V = 0.03639(15) [30] , (45) 



w = 0.84(4) [29] 



(46) 



Besides, although non-universal, let us quote e c = 
-0.867433(12) [38]. 

The results obtained from a extrapolation using only 
leading order scaling corrections were: 

• e c = -0.867397(6), u + l/u m = 1.918(26) 

(we obtained a good fit for L > L m - m = 12, with 
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X 2 /doi = 0.39/3, C.L.=94%, where "dof" stands 
for degrees of freedom and "C.L." for confidence 
level [39]). 

• £,e ,L/L = 0.5003(12), lu = 0.581(27) 
(L min = 12, x 2 /dof= 0.12/3, C.L. =99%). 

• v m = 0.714(28), u = 0.53(30) 

(L min = 8, x 2 /dof= 3.16/4, C.L. =53%). 

• r/ = 0.0391(15), u= 1.21(24) 

(Lmin = 8, x 2 /dof= 0.96/4, C.L =92%). 

The main conclusions that we draw from these fits are: 
(i) the exponents are compatible with our expectations 
from Fisher-rcnormalization, (ii) sub-leading scaling cor- 
rections are important given the tendency of the fits to 
produce a too low estimate for u> (see below) and (iii) the 
estimates from canonical exponents (obtained themselves 
by applying the high-temperature expansion to improved 
Hamiltonians [29, 30]) are more accurate than our direct 
computation in the Microcanonical ensemble. 




FIG. 2: (color online). Crossing points of the correlation 
length in lattice size units for the three dimensional Ising 
model. The error bars are in every case smaller than the point 
sizes. The values of the different quantities at the crossing as 
well as the critical exponents are shown in Table II. 

We can, instead, take an opposite point of view. If we 
take the central values in Eqs. (44,45,46), as if they were 
exact, we can obtain quite detailed information on the 
amplitudes for scaling corrections: 

• We find an excellent fit to is m (L,2L) = v m + 
Axlr* + A 2 L- 2ul , for L min = 16: x 2 /dof = 1.53/3, 
C.L =68%, with A x = 1.38(7) and A 2 = -7.6(1.1). 
This confirms our suspected strong subleading cor- 
rections. Indeed, according to these amplitudes 
Ai and A 2 , only for L ps 130 the contribution of 
the (sub-leading) quadratic term becomes a 10% of 
that of the leading one. 



• In the case of r)(L, 2L) = rj + B X L~ W + B 2 L-' 2ul , 
for L min = 8: X 2 /dof = 2.4/5, C.L =79%, we have 
B x = 0.101(10) and B 2 = 0.07(7). Subleading scal- 
ing corrections are so small that, within our errors, 
it is not clear whether B 2 = or not. 

The quite strong scaling corrections found for v m may 
cast some doubts in the extrapolation for ^L,e c /L, the 
only quantity that we cannot double-check with a canon- 
ical computation. To control this, we proceed to a 
fit including terms linear and quadratic in with 
lu = 0.84(4). We get 

^ = 0.4952(5)(7), 

with L min = 12, x 2 /3 = 2.17/3, C.L =54%. Here, the 
second error is due to the quite small uncertainty in lu. It 
is remarkable that the contribution to the error stemming 
from the error in lu is larger than the purely statistical 
one. 



VII. RESULTS FOR THE D = 2, Q =4 POTTS 
MODEL 

The Q = 4 D = 2 Potts model offers two peculiarities 
that will be explored here. First, it suffers from quite 
strong logarithmic scaling corrections. And second, it 
displays pseudo metastability [25], an ideal playground 
for a microcanonical study. 

The study of the FSS for the Q = 4, D = 2 Potts 
model [23], based on the analysis of the Renormaliza- 
tion Group (RG) equations [31], reveals the presence of 
multiplicative scaling corrections. This is one of the pos- 
sible forms that scaling corrections can take in the limit 
cij — > 0, and is a great nuisance for numerical studies. A 
very detailed theoretical input is mandatory to perform 
safely the data analysis. We shall make here an educated 
guess for the microcanonical form of the scaling correc- 
tions, based purely in ensemble-equivalence and in the 
canonical results. 

From ensemble-equivalence we expect 

e-e c ~C(L,{3 c )Af3 L , (47) 

where C(L, (3 C ) is the finite-lattice canonical specific heat 
at p c and A/? = (3^ — (3 C is the inverse-temperature 
distance to the critical point of any L-dependent feature 
(such as the temperature maximum of the specific-heat, 
etc.). We borrow from Ref. [23] the leading FSS behavior 
for these quantities: 

Thus, we have: 

e(L) - e c (oo) ~ L~ 1/2 (log L)" 3/4 . (49) 
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This result can be derived as well by considering only the 
leading terms of the first derivative of the singular part of 
free energy respect to the thermal field, </>(oc f3 — (3 C ) [23]: 



D. 



,|4/3 



(-log 



(50) 



The previous equation describes the energy of the system 
and its leading term is 



4 |0|V3 
3 ± log |0| 



but 



0«C^L- 3 / 2 (logL) 3 / 4 



(51) 



(52) 



so it is direct to obtain again Eq. (49). Hence, we are 
compelled to rephrase Eq. (27) as 



(0) L , e = L^fo (V/ 2 (logL) 3 / 4 ( e - e c )) + 



(53) 



Furthermore, from the canonical analysis [23], we ex- 
pect multiplicative logarithmic corrections to the sus- 
ceptibility (that do not get Fisher renormalized) . Fur- 
thermore, the dots in (53) stand for corrections of order 
log log L I log L and 1 / log L [23] . 

We first address in Sect. VII A the direct verification 
of Eq. (53) using the quotients method. We consider 
afterwards the pseudo-metacstability features. 



A. Scaling Plots and Critical Exponents 



onto a single curve (the deviation will be bigger for small 
L values, due to neglected scaling corrections of order 
log log L I log L and 1 / log L) [40] . A similar behavior is 
expected for the scaled susceptibility [23]: 



X 



LV^logL)- 1 /* 



(54) 



Note that £/L does not need an additional logarithmic 
factor. These expectations are confirmed in Fig. 3, spe- 
cially for the largest system sizes (that suffer lesser scal- 
ing corrections). 

We can check directly the importance of the multi- 
plicative logarithmic corrections for the susceptibility by 
comparing \ an d X as a function of Fig. 4. The 
improved scaling of x 1S apparent. We observe as well 
that the largest corrections to scaling are found at and 
below the critical point (around £/L w 1.0). 




0.4 0.6 0.: 

£(e, L)/L 



0.5 
g " 1.4 

o 

0.6 



32 
64 
128 



256 
512 
1024 



-0.5 



0.5 



, , ,1/2,, T ^3/4 

(e -e )L (logi) 



FIG. 3: (color online). Graphical demonstration of Eq. (53) 
as applied to the microcanonical D — 2, Q — 4 Potts model: 
both the correlation length in units of the lattice size (top) 
and the scaled susceptibility, x m Eq. (54) (bottom), are 
functions of the scaling variable (e — e c )Z, 1//2 (logZ,) 3 ' 4 . 

We start by a graphical demonstration of Eq.(53): £/L 
as a function of (e — e c )L 1 / 2 (log L) 3 / 4 , should collapse 



FIG. 4: (color online). Comparison of the scaling for the 
naively scaled susceptibility xL~ 7 ^ 4 (top) and for \ (bot- 
tom), as a function of the correlation length in units of the 
lattice size, for the microcanonical D = 2, Q = 4 Potts model. 

We now proceed to the numerical computation of crit- 
ical exponents. We shall use the quotients method, mod- 
ified as described in Appendix A. As it is evident from 
Fig. 5, the crossing points can be obtained with great ac- 
curacy using parabolic interpolations of the nine points 
around the estimated crossing energies, see Sect. IV. We 
checked that the results do not depend on the interpolat- 
ing polynomial degree by comparing with interpolations 
using cubic curves. We also compared with the results 
obtained using only seven points around the crossing ob- 
taining again full agreement. 

The obtained critical exponents are shown in Table III, 
we may compare them with the exact ones [26] (y = 2/3, 
a = 2/3 and 77 = 1/4): 



= 2 



V = Vn 



(55) 



Comparing with our computed exponents we obtain an 
acceptable agreement. In the case of the microcanoni- 
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cal v exponent, i/ m , after adding the correction for the 
quotients method in presence of logarithms, the agree- 
ment is fairly good. We can see a clear trend towards 
the exact result value for all the lattice sizes except for 
the biggest one (2.5 standard deviations away), which is 
probably due to a bad estimation of the huge temper- 
ature derivatives of the correlation length. In the case 
of the microcanonical r\ exponent, 7y m , which must be 
the same that the canonical one, we can see clearly the 
tendency to the analytical value 77 m = 0.25. We must re- 
mark the importance of adding the corrections described 
in Appendix A to the quotients method. 




-1.05 



-1.04 



FIG. 5: (color online). Correlation length in lattice size units 
for the two-dimensional Q = 4 Potts model. The values of 
the different quantities on the crossings for lattices L and 2L, 
as well as the corresponding estimate for critical exponents, 
are in Table III. The inset is a magnification of the critical 
region. 



B. Critical point, latent heat and surface tension 

It has been known for quite a long time that the D = 
2, Q = 4 Potts model on finite lattices show features 
typical of first-order phase transitions [25]. For instance, 
see Fig. 6, the probability distribution function for the 
internal energy, Pp(e), display two peaks at energies ed 
(the coexisting disordered phase) and e (the energy of 
the ordered phase) separated by a minimum at e*. Of 
course, since the transition is of the second order, e c is 
the common large L limit of ed, e and e*. 

We discussed in Sect. II B how the Maxwell construc- 
tion is used to estimate the canonical critical point /3 c ,l, 
as well as ed, e and the associated surface tension. This 
procedure is outlined in Fig. 7. The numerical results are 
in Table IV, where we see that (3 Cy L is a monotonically 
increasing function of L continuously approaching to the 
analytical value (3 C = log(l + y/Q) = 1.0986122 . . . [32]. 
A jackknifc method [9] is used to compute the error bars 
for all quantities in Table IV. 



20 
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-0.95 



FIG. 6: (color online). Canonical probability distribution 

function for the energy density, (e), as reconstructed from 
microcanonical simulations of the D = 2, Q = 4 Potts model 
and different system sizes. The L dependent critical point /3 c ,l 
is computed using the Maxwell rule, Sect. II B (note the equal 
height of the two peaks enforced by Maxwell's construction). 
The system displays an apparent latent heat, that becomes 
smaller for growing L, and vanish in the large L limit. 




1.09 



1.09115 



1.09105 




1.0985 



1.098 





\\ /3 c ( Exact ) " 








x 






: 



-1.06 



-1.02 



FIG. 7: (color online). Top: From the microcanonical mean 
values (/3) e ,L for the D — 2, Q = 4 Potts model, we estimate 
the size dependent canonical inverse critical temperature f3 c ,L 
(horizontal lines) for all the simulated lattice sizes, ranging 
from L — 32 (lower) to L — 1024 (upper). We show as well as 
the analytical prediction (upper horizontal line). Bottom- 
left: example of Maxwell construction for our L — 32 data. 
The e-integral of {$) e ,L — Pc.l from e to ed vanish. Bottom- 
right: zoom of upper panel showing only data for lattice sizes 
L — 256 (lower curve), L = 512 (medium curve) and L = 1024 
(upper curve). 



To perform a first check of our data, we observe that 
/3 Cj l is a typical canonical estimator of the inverse critical 
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L 


e c ,L,2L 


^,L,e c L 2L /L 








Vm 


32 


-1.04659(5) 


0.8016(5) 


1.534(6) 


1.998(10) 


0.2663(9) 


0.2334(9) 


64 


-1.04633(2) 


0.7990(3) 


1.554(8) 


1.957(12) 


0.2638(6) 


0.2360(6) 


128 


-1.04579(1) 


0.7909(3) 


1.578(5) 


1.938(7) 


0.2639(5) 


0.2398(5) 


256 


-1.04548(2) 


0.7836(5) 


1.643(12) 


1.987(17) 


0.2615(11) 


0.2402(11) 


512 


-1.04519(2) 


0.7734(9) 


1.602(31) 


1.895(42) 


0.2617(21) 


0.2427(21) 



TABLE III: Crossing points of the correlation length in lattice size units as a function of the energy for pairs of lattices (L, 2L). 
Using the original quotients method [9] we obtain the microcanonical critical exponents, shown in the columns 4 and 6, while 
the corrected ones (columns 5 and 7) are labelled with primed symbols, see Appendix A. 



L 




e G 


ed 


E x 10 5 


32 


1.0911070(20) 


-1.0175(4) 


-0.9760(2) 


0.47(2) 


64 


1.0957256(14) 


-1.0392(3) 


-0.9915(2) 


2.77(7) 


128 


1.0975150(10) 


-1.0463(3) 


-1.0062(5) 


4.10(15) 


256 


1.0981989(5) 


-1.0489(2) 


-1.0183(3) 


3.92(8) 


512 


1.0984570(3) 


-1.0490(1) 


-1.0266(2) 


3.28(11) 


1024 


1.0985539(3) 


-1.0483(3) 


-1.0325(1) 


2.09(17) 



TABLE IV: Using Maxwell construction, we compute for the 
D — 2, Q — 4 Potts model the //-dependent estimate of the 
(inverse) critical temperature /3 c ,l, the energies of the coex- 
isting ordered phase e , and disordered phase ed, as well as 
the surface tension (£)■ 



temperature. As such, it is subject to standard canonical 
FSS, where the main scaling corrections come from two 
additive logarithmic terms [23]: 



(log£) 3/4 
PcL - ft = ax L3/2 x 

log log L 

+ a 2 — ha 3 

logi 



logL 



(56) 



From our data in Table IV, we obtain a\ = —0.44(7), 
a 2 = -1.15(72), and a 3 = 2.28(26), and a good fit 
(i min = 128: x7dof = 0.28/1, C.L.=60%). 

As for the L dependence of ed and e , we try a fit that 
consider the expected scaling correction terms [23] : 



1 +a 2 



rL-V^logi)- 3 /^ 
log log L 



logL 



a 3 



logi 



(57) 



Our results for e D are: a\ Q = —2.03(20), a 2o = —1.65(27), 
and a 3o = —2.08(41), with a fair fit quality (L m in = 32: 
X 2 /dof = 2/3, C.L.=57%). On the other hand, we 
obtain for e d : aw = 2.02(14), a 2d = 0.93(37), and 
a 3 d = —2.93(34) , with a fair fit as well (L m in = 32: 
X 2 /dof = 0.84/3, C.L.=84%). These two fits are shown 
in Fig. 8. 

For the surface tension, we note in Table IV a non 
monotonic behavior. Furthermore, we lack a theoretical 
input allowing us to fit. We thus turn to a variant of the 
quotients method. Where S to follow a pure power law 




FIG. 8: (color online). System size dependent estimates for 
the energies of the "coexisting" ordered (e , blue points) and 
disordered (ed, red points) phases of the D = 2, Q = 4 Potts 
model, as a function of L" 1 ^ 2 . The lines are fits to the ex- 
pected analytical behavior Eq. (57). The horizontal line cor- 
responds to the asymptotic value, e c . 



scaling, S oc L b , exponent b would be obtained as: 



S(L 2 ) 



log^Lr)/^) 
log(Li/L 2 ) 



(58) 



The effective exponent b obtained from our data is dis- 
played in Table V. We observe that it is clearly negative 
(as it should since £ vanishes for a second order phase 
transition). An asymptotic estimate, however, seems to 
require the simulation of larger systems. 



(Li,L 2 ) 



M^) 



(32,64) 
(64,128) 
(128,256) 
(256,512) 
(512,1024) 



2.56(7) 
0.56(6) 
-0.065(60) 
-0.257(57) 
-0.650(127) 



TABLE V: Effective exponent obtained using Eq. (58) for the 
surface tension. 
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L 


Heo)/L 




X(e ) 


X(e d ) 


^canonical j 


^-canonical 
A 


32 


0.637(2) 


0.453(1) 


0.990(3) 


0.907(2) 


0.647(1) 


1.287(3) 


64 


0.732(3) 


0.396(1) 


0.995(2) 


1.025(3) 


0.545(2) 


1.310(2) 


128 


0.799(5) 


0.357(4) 


1.001(3) 


1.106(5) 


0.472(7) 


1.331(3) 


256 


0.866(6) 


0.335(3) 


1.001(5) 


1.182(6) 


0.429(5) 


1.343(5) 


512 


0.915(4) 


0.315(2) 


1.014(8) 


1.238(4) 


0.392(4) 


1.366(8) 


1024 


0.953(15) 


0.302(2) 


0.997(21) 


1.279(13) 


0.367(3) 


1.353(22) 



TABLE VI: Correlation length in units of the lattice size and the RG invariant \ defined in Eq. (54), for several L values, as 
computed in the microcanonical D — 2, Q = 4 Potts model. The chosen values of the energy density correspond to the ordered 
(e ) and disordered (e<i) phases. For comparison we also display the canonical results at /3 C obtained in Ref. [23]. 



We have just seen that, up to scaling corrections, e d L ' 
and correspond to (different) L-independent values 
of the argument of the scaling function in Eq. (53). 
Hence we expect that £(ed)/L and £(e D )/L, see table VI, 
approach non- vanishing, different values in the large L 
limit. The finite-size scaling corrections are expected to 
be additive logarithms [23] 



i 
L 



logL 



The results are 



& ) 



1.28(1) 



2.28(5) 
logL 



(L E 



32: x 2 /dof= 4.2/3, C.L.=22%), and 

= 0.159(4) -« 
L logL 



(59) 



(60) 



(61) 



(L min = 32: x 2 /dof= 3.3/3, C.L.=37%). 

A very similar analysis can be performed for the scaled 
susceptibility, Eq. (42), at ed and e . In order to deal 
with the multiplicative logarithms of the susceptibility 
we rather used % defined in Eq. (54). 

Fitting our data set to the logarithmic form 



X = A + B 



log log L 
logL 



(62) 



obtained in Ref. [23] , we obtain a good fit in the ordered 
phase energy, e : 

X(e„) = 2.41(5) -4.00(15)^|^, (63) 

logL 

(L min = 128: x 2 /dof= 3.10/2, C.L.=21%). On the other 
hand the extrapolation for the susceptibility defined in 
the disordered phase energy, e^, is a nonsensical negative 
value. 

We can also fit the data to the logarithmic form also 
used in Ref. [23]: 



X = A 



B 
logL 



(64) 



finding: 



X(e ) = 1.643(5) 



2.55(2) 
logL 



(L min = 32: x 2 = 7.44/4, C.L.=11%), and 



X(e d ) = 0.094(7) 



1.87(37) 
logL 



(65) 



(66) 



(L min = 64: x 2 /dof= 2.94/3, C.L.=37%) . For compari- 
son, we recall that Ref. [23] reports two different fits for 
X, depending of the logarithmic corrections they used: 

-canonical = 1 573(33) _ 1 .056(98) ^ L , (67) 



logL 



-canonical = 1 _ 4 54( 13 ) 



0.600(55) 
logL 



(68) 



VIII. CONCLUSIONS 

We have formulated the Finite Size Scaling Ansatz for 
microcanonical systems in terms of quantities accessible 
in a finite lattice. This form allows to extend the Phe- 
nomenological Renormalization approach (the so called 
quotients method) to the microcanonical framework. 

Our Ansatz has been subjected to a strong numerical 
testing. We have performed extensive microcanonical nu- 
merical simulations in two archetypical systems in Statis- 
tical Mechanics: the three dimensional Ising model and 
the two-dimensional four states Potts model. The two 
models present a power-law singularity in their canonical 
specific heat, implying non-trivial Fisher renormalization 
when going to the microcanonical ensemble. A micro- 
canonical cluster method works for both models, hence 
allowing us study very large system sizes (L = 128 in 
D = 3 and L = 1024 in D = 2). 

In the case of the Ising model, we have obtained precise 
determinations of the critical exponents, that, we feel, 
provide strong evidence for our extended microcanonical 
FSS Ansatz. 

For the Potts model, very strong logarithmic correc- 
tions (both multiplicative and additive) plague our data. 
Fortunately, we have a relatively strong command on 
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these corrections from canonical studies [23]. Our data 
can be fully rationalized using the scaling corrections sug- 
gested by the theoretical analysis. 
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APPENDIX A: THE QUOTIENTS METHOD IN 
PRESENCE OF MULTIPLICATIVE 
LOGARITHMIC CORRECTIONS 

The quotients method [9, 22], as been widely used in 
the past for the computation of critical exponents. Yet, 
its convergence to the large L limit is extremely slow 
in presence of multiplicative logarithmic scaling correc- 
tions. Fortunately, let us show how we can speed up 
convergence if we have enough analytical information at 
our disposal. 

Let us consider an observable O such that its FSS be- 
havior is given by (z can be either the reduced tempera- 
ture t or e — e c ) 



0(L,z) = L x °/ u (\og Vf° 



, (Al) 



then the critical exponent calculated using Eq. (31) must 
be corrected following: 



xp 
v 



xo 



log(L 2 /£i 



log 



logi] 



(A2) 



Specifically for the two dimensional four states Potts 
model the values of the logarithmic correction exponents 



are analytically known [23, 31] thus we can calculate ac- 
curately the corrections in this case. In addition, the 
susceptibility behaves as 



X 



LV'QogL)- 1 '* 



so we easily get 
rj = r\ 



1 



log 



logi; 



(A3) 



(A4) 



81og(L 2 /ii) Vlogii 
For the correlation length it is known that 

£ Hip 2 / 3 (- logi) 1/2 ; t~L-^(lo g L)^ 

(A5) 

and therefore his temperature derivative scales as 

a^^i 5/2 (logi)- 3/4 (A6) 
resulting in a v canonical exponent correction of 



1 



log 



log £5 



(A7) 



41og(L 2 /Xi) VlogLi, _ 
While for the microcanonical v exponent, u m , we use 
that 



and 



Hence, 



e ~ L~ 1/2 (logL) _3/4 
cU~L 3 / 2 (logL) 3 / 4 . 



log 



( log 



41og(L 2 /ii) °Vl°g i i 



(A8) 



(A9) 



(A10) 
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Indeed, it suffices to notice that (C(T) is the canonical 
specific heat, C oc |£| -a ), 
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Note that this microcanonical ensemble exactly matches 
the conditions in the original Fisher work [10]: the mo- 
menta are some hidden degrees of freedom in thermal 
equilibrium with the spins, and a global constraint is im- 
posed. It is also amusing to rederive the results in Sect. II 
considering r momenta per spin (in this work F — 1, 
while Lustig [5] always considered r = 3). If one takes 
the limit r — » oo, at fixed N, the canonical probability 
is recovered for the spins. 

Note that the microcanonical weight (9) is not analytical 
at each energy level of the spin Hamiltonian. 
Note that, Eq. (27) tells us that, if the energy histogram 
is double-peaked, see Sect. II B, the histogram maxima 
will tend to e c only as L^ 1 ^" 1 . 

For the 3D Ising model at criticality, ujf™ 8 = 
-0.990627(24) [27], and #[ sing = 0.2216546(2) [28], wc 
obtain for our Potts representation of the Ising model 
ec = (^ sing -D)/2 + l/(4^ sing ). 

The confidence level is the probability that \ 2 would be 
bigger than the observed value, supposing that the sta- 
tistical model is correct. As a rule, we consider a fit not 
good-enough whenever C.L.< 10%. 

We obtain the exact e c in the thermodynamic limit from 
P c = log(l + y/Q) [32], and u c = -(1 + Q~ 1/2 ) [26] by 
applying e c = u c + l/(2/3 c ) 



